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We review recent attempts to address the cosmological constant problem and the late-time accel- 
eration of the Universe based on braneworld models. In braneworld models, the way in which the 
vacuum energy gravitates in the 4D spacetime is radically different from conventional 4D physics. 
It is possible that the vacuum energy on a brane does not curve the 4D spacetime and only affects 
the geometry of the extra-dimensions, offering a solution to the cosmological constant problem. We 
review the idea of supersymmetric large extra dimensions that could achieve this and also provide 
a natural candidate for a quintessence field. We also review the attempts to explain the late-time 
accelerated expansion of the universe from the large-distance modification of gravity based on the 
^ ' braneworld. We use the Dvali-Gabadadze-Porrati model to demonstrate how one can distinguish 

^s^j , this model from dark energy models in 4D general relativity. Theoretical difficulties in this approach 

are also addressed. 



I. INTRODUCTION 

The cosmological constant problem is a long-standing problem in physics Particle physics predicts 
the existence of the vacuum energy density which is related to the fundamental scale of the theory, like the 
electroweak scale, p vac ~ (TeV) 4 . This is typically more than 50 orders of magnitude larger than the observed 
value, pa ~ (10" 3 eV) 4 . Before the discovery of the accelerated expansion of the Universe, physicists tried to 
answer this question by seeking a theory that predicts the cosmological constant should be zero. However, the 
discovery of the accelerated expansion of the Universe makes this answer insufficient 0, [H, [f| . Now, we should 
explain why it is non-zero and yet it is so small. Moreover, there is a coincidence problem. The cosmological 
' constant dominates the energy density of the Universe only recently. If the cosmological constant is really a 

constant, we should explain why now, does it become dominant. 

One direction to answer these questions is to appeal to the anthropic principle If the cosmological 
constant is too large, the accelerated expansion started too early and it prevents structure from growing 
and we cannot exist. On the other hand, a universe with negative comsological constant re-collapses. Then 
observers will only exist within a tiny anthropic range of cosmological constant (see for example @). This 
idea is strengthened by the discovery in string theory that there are millions of low-energy vacua in the 
theory (the string theory landscape) Q- It is argued that we might need the anthropic principle to select 
the low-energy vacuum. However, many theorists still hope to explain the problem without invoking the 
existence of ourselves in the Universe. Although significant efforts have been devoted to this attempt, we 
still have not succeeded yet to provide convincing models. However, the rapid progress of string theory 
has provided a new perspective for solving this problem. In this review, we focus on the attempts of using 
higher-dimensional gravity and branes to address the problem. 

String theory is formulated in a 10D spacetime. On the other hand, our observed Universe is a 4D 
spacetime. Thus there should be a mechanism to hide the extra dimensions. The conventional idea is to 
compactify the extra dimensions by the Kaluza-Klcin (KK) mechanism. The size of the extra dimensions L 
should be small, L < TeV -1 , in order not to spoil the success of the standard model of particle physics that 
is formulated in a 4D spacetime. Below the energy scale determined by the size of the extra dimensions, L~ , 
the universe looks completely 4D if the radius of the extra dimensions is stabilized. Recently, a completely 
new way of hiding the extra dimensions has been proposed. This is the brane world mechanism where 
matter fields are confined to a 4D membrane in a higher dimensional spacetime (see [8j] for a review). Only 
gravity and non-standard model particles can propagate into the whole higher-dimensional bulk spacetime. 
In this picture, the size of the extra dimensions can be much larger than that in the conventional KK 
compactification. In fact, the size of the extra dimensions could even be infinite. If the bulk is a spacetime 
with a negative cosmological constant, that is, an Anti-de Sitter (AdS) spacetime, it is shown that gravity 
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behaves like 4D on scales larger than the AdS curvature length, even if the size of the extra dimensions is 
infinite 9]. Another way is to introduce induced gravity on a brane [l(| (see [ll[ for an early attempt). If 
we assume there is an Einstein-Hilbert term on a brane, 4D gravity is recovered, in this case, on small scales 
even if the bulk is an infinite Minkowski spacetime. In these braneworld models, the behaviour of gravity can 
be dramatically different from the 4D theory, providing a new perspective to solve the cosmological constant 
and the dark energy problem. 

This article will review several approaches to address the cosmological constant and the late-time accelera- 
tion problem based on braneworld gravity. Firstly, we explain the attempts to address the 'old' cosmological 
constant problem — why the vacuum energy is incredibly small compared with the prediction of particle 
physics. These attempts exploit the modification of 4D gravity in the braneworld and change the way in 
which the vacuum energy gravitates in a 4D spacetime. Secondly, we introduce an idea to explain the late- 
time acceleration without introducing the cosmological constant. This idea also relies on the modification of 
gravity on large scales based on the braneworld idea. 

In section II, we give a brief introduction to braneworlds. In section III, the attempts to solve the old 
cosmological constant problem are discussed. In section IV, the idea to realize late-time acceleration without 
introducing a cosmological constant is explained. Section V is devoted to conclusions. 

II. BRANEWORLD MODELS 

The idea that ordinary matter fields are confined to a lower-dimensional domain wall was proposed in 
the 1980's [HI, EH- It was shown that fermion fields can be confined to a field theoretic domain wall. The 
progress in string theory, especially the discovery of D-branes, has revived these attempts [l4j . The D-brane 
is defined by a membrane on which end-points of open strings lie. At the end-points of open strings, gauge 
fields can be attached. Then gauge fields are confined to the D-brane. On the other hand, closed strings that 
contain the graviton can propagate into the whole bulk. Then there arises a braneworld picture where usual 
matter fields are confined to a brane while gravity propagates into the whole bulk spacetime. A schematic 
picture of the braneworld is shown in Fig. 1. Based on this idea, several simplified braneworld models have 
been proposed that capture the basic features of the braneworld, yet in which we can address many important 
problems from a new perspective. 




FIG. 1: A schematic picture of the braneworld. From [1]. 
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FIG. 2: Constraints on Yukawa violations of the gravitational 1/r potential , V(r) oc (l/r)(l + aexp(— r/A)). The 
shaded region is excluded at the 95% confidence level. From [f8l ]. 



A. Arkani-Hamed-Dimopoulos-Dvali model 

An interesting possibility in braneworld models is that some of the extra dimensions can be large [l5|, [T5 , 
[13]. In a conventional picture, extra dimensions are rolled up small so that we never observe them. More 
precisely, in order not to spoil the success of the standard model of particle physics that is formulated in a 
4D spacetime, the size of the extra dimensions should be smaller than TeV -1 ~ 10~ 19 m. However, in the 
braneworld, the standard model particles are confined to the 4D brane. Thus we do not need to worry about 
this constraint. The gravitational interactions are very weak and the 4D behaviour of the Newtonian force 
is only verified down to 44/itm [18j . Thus the size of the extra-dimensions is allowed to be as large as 44/im. 

This opens up a new perspective to solve another serious problem in particle physics, namely the hierarchy 
problem: why the gravitational interaction is so weak compared with the other interactions. The answer 
could be that the gravitational field of an object on a brane leaks out into the large extra-dimensions and 
this leakage weakens the gravitational interactions on a brane. The gravitational potential generated by an 
object with mass M is given by 

*W = <r < l), ( 2 ) 

where L is the size of the (D — 4) dimensional extra-dimensions. Then the 4D gravitational constant is given 
in terms of the higher-dimensional gravitational constant as 

G 4 = ^£_ M 4 = M(?- 2 >/W- 4 )/ 2 , (3) 

where 8nG D = M D [D ~ 2) . Then even if the fundamental scale of gravity Mo is TeV, the 4D gravitational 
constant can be 10 19 GeV as long as L is appropriately large. For example, for D = 6, the current constraint 
on the deviation from the gravitational inverse-square law L < 44/im implies Me > 3.2 TeV. 
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Randall-Sundrum model 



The most difficult problem in the braneworld in terms of gravity is the inclusion of the self-gravity of the 
branes. In the ADD model, the self-gravity of the branes is implicitly neglected. The model proposed by 
Randall and Sundrum (RS) offers a consistent framework to deal with higher-dimensional gravity including 
the self-gravity of the branes 0] . They consider a 5D spacetime described by the action 



S = 2k| ' d5x "/-9( (5)R - 2A ) _ a I d i x^f^+ I d i xC mi (4) 



where k 2 = 8ttG^ and C m represents the matter lagrangian confined to a brane. The introduction of the 
singular objects enforces the junction condition (Israel junction condition) at the location of the brane. The 
junction condition relates the extrinsic curvature at the brane to the energy momentum tensor localized on 
a brane. By solving the 5D bulk spacetime and imposing the junction condition at the brane, the solution 
for the gravitational field on the brane is obtained. The simplest solution is a solution with a Minkowski 
brane. The 5D metric is given by 

ds 2 = dy 2 + exp(-2|y|/£)^ ! ,cfe'W. (5) 

A brane is located at y = and the reflection symmetry (Zi symmetry) across the brane is imposed. 
The exponential 'warp factor' is an essential ingredient of the model. Even if the physical size of the fifth 
dimension is infinite, low-momentum gravity is confined near the brane due to the curvature of the bulk 
spacetime and 4D grav ity is recovered. It is shown that the solutions for weak gravity at large distances 
r ^> i are given by [19| 

ds 2 = -(1 + 2^)dt 2 + (1 + ^)S i3 dx l dx j , (6) 

2G 4 M ( ^ , 2£ 2 \ (J) _ 2G A M f I 2 
3r 2 / ' r \ ir' A 



where k 2 = 87rG4 and is determined by 

G 4 = G 5 l (8) 

Comparing this with Eq. |(5J), we notice that £ acts as the effective size of the extra dimension. Thus the RS 
model provides an 'alternative to compactification'. 

Despite the remarkably simple setup of the model, gravity in this model is incredibly complicated. For- 
tunately, for a homogeneous and isotropic brane, the generalized Birkoff theorem ensures that the bulk 
spacetime is AdS spacetime or AdS-Schwarzchild spacetime. Then the Friedmann equation on the brane is 
easily derived as [lo, HH, [22j 

H 2 = ^± 4 4 2 C 

3 3 H W a 4 ' y ' 

where 

As 4£ 4a 

2 12 ' 4 6 K ' 

The constant C is proportional to the black hole mass in the bulk. In accord with weak gravity, cosmology 
also shows the transition from 4D to 5D. At high energies HI > 1 where the horizon size H~ l is smaller 
than the effective size of the extra-dimension i, the Friedmann equation is significantly modified and H oc p. 
At low energies Hi < 1, we recover the 4D Friedmann equation. 

III. COSMOLOGICAL CONSTANT PROBLEM IN THE BRANEWORLD 
A. Self-tuning 5D braneworld 



The relation between the vacuum energy and the effective cosmological constant on a brane is different 
from that in the usual 4D theory. In the RS braneworld, the vacuum energy in the brane a is not directly 
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related to the cosmological constant A4 on the brane in the effective Einstein equation as in Eq. (jTUJ) . In the 
RS braneworld, there should be a cancellation between the 4D and 5D contribution of the vacuum energy in 
order to have a vanishing cosmological constant on the brane. This requires a fine-tuning for the parameters 
in the action. Instead of having the cosmological constant in the bulk and tension on the brane, let us 
consider a scalar field with potentials 0, [24| • The action is given by 



1_ f , K , ( _ 4 



d b xV=g~ ( R - ^{d^f - V{cj>) ) - I d 4 xV^/(0). ( 1 1 ) 



The potentials can be taken as 



V(<j>) = A n exp(a^), f(</>) = V exp(ty). (12) 



With this choice, the action describes a family of theories parametrized by Vq, Ao,a and b. For simplicity, 
we take Ao = 0. We look for a solution with a Minkowski spacetime on a brane. The 5D metric is given by 

ds 2 = dy 2 + e 2A ^ v) r] liv dx l, dx fi . (13) 

The 5D Einstein equation gives the relation between the warp factor A(y) and the scalar field <fi(y) 

<P>(y) = ±~A'(y). (14) 
The solution for <f> in the bulk is then obtained as 

4>(y) = ^log \^M 5 y + cA +d 1} y < 0, (15) 
<t>{y) = — log (^M 5 y + c 2 ^ +d 2 , y>0, (16) 

where Cx,c%,di and d 2 are integration constants. The continuity of <f> determines d 2 - Then the junction 
conditions for the scalar field and the warp factor determine c\ and ci in terms of 6, Vq and d\ if b =/= ±4/3. 
This means that for a scalar coupling given by b, there is a Minkowski solution on a 4D brane for any value 
of the brane tension Vq. This is the idea of the 'self-tuning'. The vacuum energy in a 4D brane is cancelled 
by the integration constants in the solutions, not by the parameters in the original action. Thus this is not a 
fine-tuning. The hope is that the solution in the bulk adjusts itself so that the contribution from the vacuum 
energy on the brane is exactly cancelled. 

Although the idea of self-tuning is very attractive, there are several problems in the original proposal 
[25l [26l [27l [28| . Firstly, there is a naked singularity in the above model with a scalar field. Any procedure 
that regularizes the singularity in the solutions would cause the re-introduction of the fine tuning. There 
is also a problem of stability. In the case of vanishing potential in the bulk, the static solution is unstable, 
leading to a singularity. A modified version of the model using the bulk black hole to hide the singularity 
inside the horizon was proposed [2^ |. but it was argued that this model also cannot avoid the fine tuning 

M. 



B. 6D braneworld 



Another approach to realize the self-tuning is to consider a 6D bulk spacetime [31]. The action is given by 

1 



/ '!"■>■ v .</ ( - A 6 - -F ab F ab ) . (17) 

where the gauge field F a b is required to stabilize the size of the extra dimensions. We decompose the 
coordinates into four macroscopic dimensions and the two extra dimensions. The metric is taken as 

ds 2 = r\ liv dx ll dx v + j lj dx l dx j . (18) 





FIG. 3: Removing a wedge from a sphere and identifying opposite sides to obtain a football geometry. Two equal- 
tension branes with conical deficit angles are located at either pole; outside the branes there is constant spherical 
curvature. From [311 ] . 



The gauge field is taken to consist of magnetic flux threading the extra dimensional space so that the field 
strength takes the form 



Ft, 



(19) 



where Bq is a constant, 7 is the determinant of 7^ and is the antisymmetric tensor normalized as 
ei2 = 1. All other components of F a b vanish. A static and stable solution is obtained by choosing the 
extra-dimensional space to be a two-sphere 



r y i jdx l dx J = al(d9 2 + sin 2 6d(p 2 ). 
The magnetic field strength Bq and the radius ao are fixed by the cosmological constant 



B 2 = 2A 6 , a; 



Ml 

2A R 



(20) 



(21) 



It should be noted that Bq has to be tuned so that a Minkowski spacetime is induced in 4D. Now we add 
branes to this solution. The brane action is given by 



a d sc-v/— ' 



The solution for the extra dimensions is now given by 

jijdx'dx 3 = al(dd 2 + a 2 sin 2 9dip 2 ), 



where 



27rM| ' 



Ml 

2A fi ' 



(22) 



(23) 



(24) 



a) in 



The coordinate ip ranges from to 27r. Thus the effect of the brane makes a deficit angle S = 2w(l 
the bulk. This is a 6D realization of the ADD model including the self-gravity of branes. 

The most interesting feature of this solution is that the 4D geometry is independent of the brane tension 
a. The tension enters only in the deficit angle and not the radius ao and the magnetic field Bq that need 
to be tuned to obtain a flat 4D spacetime. Thus the vacuum energy on the brane does not gravitate in the 
4D spacetime but merely changes the geometry of the extra dimensions. Thus the outcome is similar to the 
self-tuning solutions discussed in section III. A. It should be noted that the cosmological constant problem is 
not fully solved even if this idea works. In order to obtain a flat spacetime, we need to tune the magnetic field 
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and the bulk cosmological constant as in Eq. (|2ip . However one could hope that some kinds of symmetry 
like supersymmetry in the bulk can ensure this tuning. 

However, there have been objections to the self-tuning in this model 32]. Consider that a phase transition 
occurs and the tension of the brane changes from o\ to 02- Accordingly a changes from ai = 1 — (7i/(27rM 4 ) 
to ot2 = 1 — oil (2irM 4 ). The magnetic flux is conserved as the gauge field strength is a closed form, dF = 0. 
Then the magnetic flux which is obtained by integrating the field strength over the extra dimensions should 
be conserved 



$s = 47raiB ,i = 47ra 2 -B 0i 2- (25) 

The fine-tuning of A$ and Bq, Eq. (|2"Tj) , that ensures the existence of Minkowski branes cannot be imposed 
both for B — £? ,i an d B — £>o,2 when a\ 7^ This becomes clear if we rewrite the conditions Eq. ([2"T]) 

as 

2 _ (®b\ 2 A 6 



<* = hr5 lir (26) 



4tt J M 6 4 



The left-hand side changes by the phase transition but the right-hand side cannot change. Moreover, the 
quantization condition must be imposed on the flux <£> b ■ Then if the condition Eq. (|21[) is satisfied for some 
value of a, it will not be satisfied by neighbouring values. Thus after the phase transition, the 4D spacctime 
cannot be static 1331] - 



C. Supersymmetric large extra dimensions 



In the Einstein-Maxwell theory discussed in section III.B, the tuning between the magnetic flux and the 
cosmological constant in the 6D spacetime, Eq. (f2"Tj) . was necessary to obtain the flat 4D spacetime. This 
was the origin of the difficulty in realizing the self-tu ning . To evade this problem, the Supersymmetric Large 
Extra Dimensions (SLED) model was proposed (see |34l. |35| for a review). This is a supersymmetric version 
of the 6D model and the action is given by [36| 

f 



S 



2Kg 



R-d 



5M, 



MN 



(27) 



There exists a solution where the dilaton </> is constant, (j> — 
theory is a solution just by replacing Ag — > Age^ and B% 
by the condition that the potential for 4> has minimum 



0, and the solution in the Maxwell-Einstein 
The constant value <fio is determined 



B%e 



1 9 
-Ble 

2 



V'(<h) = --B'te-* + Aee^ = n. 
This is exactly the condition to have a flat geometry on the brane (see Eq. (|2"Tj) ) 

Ble'* = 2A 6 e*°. 



(28) 



(29) 



Thus unlike the Einstein-Maxwell system, one might not need a tuning condition in the bulk. In fact, The 
known solutions in this model which have maximally symmetric 4D metric all have vanishing vacuum energy. 

Again there were several objections to this version of the self-tuning [32|, [3?|]- It is possible to derive the 
4D effective theory by putting the metric in the form 



ds 2 = g^{x)dx^dx" + M^ 2 e - 2 ^ x > (dr 2 + sin 2 rd6 2 ), 
and assuming cj> — cf>(x). The potential which results from the two scalar fields is [32j 

2A 6 , 



V(iPA)=M- 4 e°*U(a 1 ) 



2M 2 e- ai 



(30) 



(31) 



where <j\ = 2%j) + <fi and a% — 2ip ~ (p. Unlike the Einstein-Maxwell theory, &i ensures that U{<j\) vanishes at 
the minimum of the potential. We should note that 02 is related to the classical scaling property of the model. 



8 



The 6D equation of motion is invariant under the constant rescaling gmN —> e u 9mn an d e * ~~ > and the 
lagrangian is scaled as C — > e 2u> C. The modulus 02 can be identified as the one associated with this scaling 
property. Thus the flatness of the 4D spacetime is ensured by the scaling property of the theory. However, 
this eventually leads to the same tuning condition (|29[) as in the Einstein-Maxwell theory. Then we can 
apply the same arguments as in the previous section. Suppose that the tension of the brane changes. Flux 
conservation (and flux quantization) means that the tuning condition cannot be maintained and U(ai) =/= 0. 
What happens would be that 02 acquires a runaway potential and the 4D spacetime becomes non-static. 

A caveat in this argument is that the metric ansatz (|30|) is restrictive. In fact, there is a class of static 
solutions where there is a warping in the bulk. The solution has the form (38l . l39j 

dsl = W{r)Yn^dx»dx v + al{W{rifdr? + dd 2 ), <f> = </> + ilnW(rj) + 2A 3 ry, (32) 

where W{rf) is the warp factor. If both branes, at the north pole and the south pole, have the same tension, 
the warp factor becomes trivial. However, if the tensions are not equal, there is a warping. For A3 ^ 0, the 
metric near the branes no longer corresponds to that of a simple conical singularity. These solutions cannot 
be described by the ansatz ([3"0| . Thus one can still hope that the solutions will go to these solutions after 
a change of tension. An unambiguous way to investigate this problem is to study the dynamical solutions 
directly in the 6D spacetime. However, once we consider the case where the tension becomes time dependent, 
we encounter a difficulty to deal with the branes. This is because for co-dimension 2 branes, we encounter 
a divergence of metric near the brane if we put matter other than tension on a brane. Hence, without 
specifying how we regularize the branes, we cannot address the question what will happen if we change the 
tension. Is the self-tuning mechanism at work and does it lead to a 4D static solution? Or do we get a 
dynamical solution driven by the runway behaviour of the moduli field? There was a negative conclusion 
on the self-tuning in this supersymmetric model for a particular kind of regularization [37] . However, the 
answer could depend on the regularization of branes and the jury remains out. It is important to study the 
time-dependent dynamics in the 6D spacetime and the regularization of the branes [10, SI SI) EH • 

If the self-tuning mechanism works, then we should seek an explanation for the accelerated expansion today. 
The supersymmetry in the bulk would also provide a very interesting mechanism (see [34l l35j and references 
therein for detailed discussions). Supersymmetry is supposed to be broken at least at the electroweak scale 
M w . Then in the 4D spacetime, this gives a vacuum energy of the order p ~ as the cancellation 
between the contribution to the vacuum energy from boson fields and fermion fields ceases to exist at M w . 
However, if the self-tuning mechanism is at work, this vacuum energy does not give any contribution to the 
cosmological constant on the brane. However, the breakdown of supersymmetry is mediated to the bulk at 
least gravitationally. Then there arises a supersymmetry breaking scale in the bulk given by 

Interestingly, this scale is related to the size of the extra dimensions. If we want to solve the hierarchy 
problem between the Planck scale and the electroweak scale, M$ should be of the order M w . Then from the 
relation between Mq, M4 and the size of the extra-dimension L, Eq. ([3]), the supersymmetry breaking scale 
in the bulk is given by 

M sb = I (34) 

If L is 10/im, we get the correct order of magnitude for the cosmological constant if p\ ~ Mf b , In order to 
confirm this expectation, we should compute the effective potential for the radion which describes the size 
of the extra dimensions generated by supersymmetry breaking. The potential for the radion obtained by 
integrating out the bulk loops is given by 

V(L) = + ||(log(M 6 2 i 2 ) + C). (35) 

The calculation of C2 depends on the details of the spectrum of the theory at M w and C2 = is critical for 
this model to work. If C2 = 0, the potential leads to a natural realization of the quintessence model where 
the radion L acts as a quintessence field. 

Thus SLED gives a consistent framework to address the cosmological constant problem and the dark 
energy model provided that the self-tuning mechanism works and the supersymmetry breaking on the brane 
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generates the desired potential for the radion, V(L) ~ L~ 4 . In SLED, the 6D Planck scale is supposed 
to be M w and the size of the extra dimensions today are L ~ lOum. This leads to a lot of interesting 
phenomenology in local tests of gravity, collider physics and so on [34], [H| . 

IV. LATE-TIME ACCELERATION IN THE BRANEWORLD 

A new twist to the cosmological constant problem is the late time acceleration of the Universe. The simplest 
way to realize this is to assume that a tiny amount of the cosmological constant is left after cancelling the 
vacuum energies. But the vacuum energy is typically more than 50 orders of magnitude larger than the 
observed value of the cosmological constant. Thus this is an incredible fine-tuning. Moreover, if the self- 
tuning idea works and the vacuum energy does not gravitate, it is in general difficult to realize the accelerated 
expansion of universe (see however the SLED proposal discussed in section III.C). Alternatively, it is possible 
that there is no cosmological constant but that large-distance modification of GR accounts for the late-time 
acceleration. The braneword gravity provides a natural framework for the study of this possibility. For 
example in the model proposed by Dvali, Gabadadze and Porrtati (DGP), 4D GR is modified on large scales 



1C | . It is in fact possible to realize the accelerated expansion of the universe without a cosmological constant 
44, EH]. This solution is known as the self-accelerating universe. We should note that in these attempts, 
we do not solve the old cosmological constant problem. In addition, in the DGP model, the coincidence 
problem is not solved and we should introduce a fine-tuned dimensional parameter related to the scale of 
the cosmological constant, pa ~ 10~ 3 eV. However this is a novel alternative to dark energy models in GR 
and it gives a new perspective to approach the problem. 

A. Dvali-Gabadadze-Porrati model 

In the DGP model, gravity leaks off the 4D Minkowski brane into the 5D bulk Minkowski spacetime at 
large scales. The 5D action describing the DGP model is given by 

1 

c 5 



S= ^ / d b x^R+ — / dW^T {4, R- I d A x^^C m . (36) 



Instead of having the bulk cosmological constant and the tension on a brane as in the RS model, there is an 
induced Einstein-Hilbert term on the brane. 

On small scales, gravity is effectively bound to the brane and 4D Newtonian dynamics is recovered to a 
good approximation. The transition from 4D to 5D behaviour is governed by a crossover scale 



2k 



2 ' 



(37) 



The weak-field gravitational potential behaves as 



*H T ~\ ? r r<rc ' (38) 



Unlike the RS model, gravity becomes 5D at large distances. The DGP model was generalized by Deffayet 
to a Friedman-Robertson- Walker brane in a Minkowski bulk (45| . The energy conservation equation remains 
the same as in general relativity, but the Friedman equation is modified: 

p + 3H(p + p) =0, (39) 

-=^ 2 -^V (40) 

T c 3 

The modified Friedmann equation shows that at late times in a CDM universe with p oc a~ 3 — - > 0, we have 

H — > = — ■ (41) 
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Since Hq > -Hoc, in order to achieve acceleration at late times, we require r c > Hq 1 , and this is confirmed 
by fitting SN observations [5l|. Like the LCDM model, the DGP model has simple background dynamics, 
with a single parameter r c to control the late-time acceleration. 

On small scales, the Newtonian potential behaves as 4D. The Friedmann equation also shows that the 
universe behaves as 4D at early times, Hr c ^> 1. However, the recovery of GR is very subtle in this model 
[46| . In fact, although the weak- field gravitational potential behaves as 4D on scales smaller than r c , the 
linearized gravity is not described by GR. This is because there is no normalized zero-mode in this model 
and 4D gravity is recovered as a resonance of the massive KK gravitons. The massive graviton contains 5 
degrees of freedom compared with 2 degrees of freedom in a massless graviton. One of them is a helicity- 
polarization. Due to this scalar degree of freedom, linearized gravity is described by Brans-Dicke (BD) 
gravity with vanishing BD parameter in the case of Minkowski spacetime. Thus this model would be excluded 
by solar system experiments. Howe ver, the non-linear interactions of the scalar mode becomes important 
on larger scales than expected 0, H3, [H, H§|. Let us consider a static source with mass M. Gravity 
becomes non-linear near the Schwarzshild radius r g = 2GM . However, the scalar mode becomes non-linear 
at r* = (r^r;?) 1 / 3 (the Vainstein radius) which is much larger than r g if r c ~ Hq . In fact, for the Sun r* 
is much larger than the size of the solar system. A remarkable finding is at once the scalar mode becomes 
non-linear, GR is recovered. This non-linear shielding of the scalar mode is crucial to escape from the tight 
solar system constraints. Fig. [5] summarizes the behaviour of gravity in the DGP model (see [50( for a review 
on the DGP model). 

r, r r 
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FIG. 4: Summary of the behaviour of gravity in the DGP model. At large scales r > r c , the theory is 5D. On 
small scales r < r c , gravity becomes 4D but the linearized theory is described by a Brans-Dicke theory. This affects 
the large scale structure (LSS) and the Integrated Sachs-Wolfe (ISW) effect and its cross-correlation to LSS. Below 
the Vainstein radius r < r*, the theory approaches GR. This transition can be probed by weak lensing and cluster 
abundance as the non-linear dynamics is important for these measures. The solar system tests also provide constraints 
on the model in the 4D Einstein phase. 



B. Observational constraints on the self-accelerating universe 



The self-accelerating universe provides useful example where we can study how various observations can 
be combined to test the model. It also provides a possibility to find a failure of GR at cosmological scales. 
A key is the complicated behaviour of gravity. We have various cosmological observations that cover various 
scales. Then combining the various data sets, we can probe the complicated behaviour of gravity in this 
model. A central question is whether we can distinguish the DGP model from dark energy models in GR. 



1. Expansion history 



The first question is whether one can distinguish between the self-accelerating universe and the simple 
ACDM model in GR. Both models have the same number of parameters and phenomenologically both 
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theories have the same simplicity. In terms of density parameters, the Friedmann equation in the ACDM 
model is given by 



M 



n A + n K = 1. 



On the other hand in the DGP, we have [5_1| 



K 



(42) 



(43) 



where we defined 



(44) 



In order to constrain the density parameters, we can combine data from supernovae, the cosmic microwave 
background shift parameter, and possibly the baryon oscillation peak [52|, [H, [H, [H, [H, H3] • Interestingly, 
the current observations already give us a hint how we can distinguish the models. While the ACDM model 
fits the three data sets comfortably, there is some tension between the data and DGP (Fig. [5]) [53]. It is 
suggested that a slightly open universe can fit the data set better in the DGP (Fig. [5]) [541 ]. 




0.5 




0.5 



FIG. 5: Joint constraints [solid thick] on DGP models from the SNe data [solid thin], the BO measure A [dotted] 
and the CMB shift parameter S [dot-dashed]. The left plot uses SNe Gold data, the right plot uses SNLS data. The 
thick dashed line represents the flat models, Q,k = 0. From [53| . 



Note that the baryon acoustic oscillation measure requires the knowledge of the power spectrum thus the 
knowledge of perturbations. Precisely speaking, the analysis must be redone for the DGP model. We expect 
that only small corrections are involved, but this problem must be addressed. The conclusion also seems to 
depend on the data set for supernovae (Table I and II) [53| . This is also true using the latest results from 
the ESSENCE and SNLS supernova data set and the Riess 07 Gold set (Fig. [7]) [56] . 



In the future, precision data will enable us to distinguish between the DGP and the ACDM more clearly. 
Fig. [8] shows the prediction of the baryon acoustic peak oscillation observed by a future WFMOS survey 
which is assumed to contain 2.1 x 10 6 galaxies, over 200 deg 2 , at 0.5 < z < 1.3 [58j. Clearly the difference 
between the two is much larger than the error bars. 
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FIG. 6: The A\ 2 between the best fit flat and open DGP versus that of a flat ACDM model. The Gold supernova 
(SN) data set is used in the top panel and the SNLS SN data set is used in the bottom panel. The DGP model 
requires curvature and a high Hubble constant. With the addition of Key Project ( KP) direct Hubble constant 
measurements, open DGP is a marginally poorer fit to the data than flat ACDM. From [54| 





best-fit 
acceleration 
parameter 


best-fit 
density 
parameter 


best-fit 
curvature 
parameter 


x 2 

value 


DGP 
LCDM 


fire = 0.125 
Qa = 0.730 


fi m = 0.270 
Qm = 0.285 


CIk = +0.0278 
SIk = -0.0150 


185.0 
177.8 



TABLE I: Best-fit parameters from SNe-CMB shift-Baryon Oscillation constraints, and \ 2 values, for the DGP and 
LCDM models. The Gold data is used for the SNe. From [H 



2. Linear growth of structure 

Although the DGP model can be distinguished from the ACDM model, background tests will never 
distinguish the DGP model from dark energy models in GR. This is because there always exists a dark 
energy model in GR that has exactly the same expansion history as in DGP. In fact as far as the background 
evolution of the Universe is concerned, the DGP is equivalent to the dark energy model whose equation of 
state is given by [59| 

w = —i 7^ — 7~T- ( 45 ) 

l+O m (a) v ; 

For small red-shift, this is well fitted by w = wq + w a (l — a) where wq — —0.78 and w a — 0.32 if f2 m = 0.3 
today [60j. Then we cannot distinguish the DGP from the dark energy model in GR. 

However, even if the background dynamics is the same, this does not mean that the dynamics of pertur- 
bations is the same. Koyama and Maartens obtained the solutions for metric perturbations on sub-horizon 



scales by consistently solving the 5D perturbations under quasi-static approximations 6l|. Scalar metric 
perturbations are given in longitudinal gauge by 

ds 2 = -(1 + 2^)dt 2 + a 2 (l + 2<$>)5 lJ dx l dx 1 , (46) 
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best-fit 
acceleration 
parameter 


best-fit 
density 
parameter 


best-fit 
curvature 
parameter 


x 2 

value 


DGP 

LCDM 


n rc = 0.130 
Q A = 0.740 


Q m = 0.255 
Q m = 0.270 


n K = +0.0300 
Q.K = -0.0100 


128.8 
113.6 



TABLE II: As in Table I for the Legacy SNe data. From [53j 

DGP n=1 ACDM (n=2) 





FIG. 7: Comparison between the results of fitting DGP and ACDM to the SNLS and ESSENCE supernova data 
set (filled in 68%, 95% and 99% confidence regions) and the Riess 07 Gold set (dotted lines). The solid black line 
corresponds to spatially flat universes. From [56f ] . 



and the perturbed energy-momentum tensor for matter is given by 



—Sp aSq t ,, 
—a~ 1 5q' t Sp5 t j 



W = ( ~Z i " r I • (47) 



2 $ = 4ttG 4 [l-—)p5, (48) 



2 * = -4ttG 4 [ 1 + ) p5, (49) 



The solutions for the brane metric perturbations are [6l| 
where 

P = 1 - 2r c i/ ^1 + ^ J , (50) 

and 

5 = dp- 3H5q. (51) 

This agrees with the results obtained by Lue, Scoccimarro and Starkman. They find spherically symmetric 
solutions by closing the 4D equations using an anzatz for the metric and checking in retrospect that the 
obtained solutions satisfy regularity in the bulk. It was shown that the solutions (|4"5]) and (|4T)|) are uniquely 
determined by the regularity condition in the bulk within our approximations. 

The modified Poisson equation |48|) shows the suppression of growth. The rate of growth is determined 
by S, and for CDM, 

k 2 

S + 2HS = , (52) 

a z 
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k/(hMpc"') 

FIG. 8: Theoretical predictions for din P(k)/d\n k assuming a sample WFMOS, where P(k) is the power spectrum of 
galaxies. The squares with error bars are evaluated with the simple simulation of the power spectrum for the ACDM 
model. The asterisks are the DGP model, but the error bars, which are almost the same as that of the ACDM model, 
are omitted for simplicity. Theoretical curves are the DGP model (dashed red curve) and the ACDM model (solid 
black curve). The parameters are n s — 0.95, fib = 0.044, Q. m = 0.27 and the linear bias is taken as bo = 1-5. From 



which leads to 

S + 2H5 = 4ttG 4 (l+^Jp5. (53) 

Thus the growth rate receives an additional modification from the time variation of Newton's constant 
through (3. 

In Fig. [9l we show the linear growth factor S/a for the DGP model, and compare it with ACDM and with 
the GR dark energy model whose background evolution matches that of the DGP model. We also show 
the incorrect DGP result, in which the inconsistent assumption of neglecting 5D perturbations is effectively 
adopted (62J- This inconsistent assumption has been made in various treatments but it leads to unreliable 
results. The correct equations for subhorizon density perturbations are crucial for meaningful tests of DGP 
predictions against structure formation observations. This highlights the fact that the growth rate is very 
sensitive to the modification of gravity. 

There are several observations that can probe the growth of structure. Weak lensing measures the deflection 
of light generated by matter fluctuations (see [63| for a review). The deflection potential is given by 



We can relate <f> to the matter overdensity 5: 



$ + (54) 



inG4a2 -S. (55) 



k 2 



Interestingly, this formula in the DGP is the same as the one in GR. However, the change of the growth 
rate leads to a different prediction of weak lensing. We should note that current observations measure weak 
lensing sourced by matter fluctuations in the non-linear regime. The solutions (|48p can be applied only to 
linear perturbations and there is no justification to use the linear growth rate and predict the non-linear 
power spectrum using the mapping formula developed in GR. We will come back to this issue in section 
IV.B. 

Another probe is the Integrated Sachs- Wolfe (ISW) effect. This is determined by the time variation of 
the deflection potential cb. On large scales, we should deal with the truly 5D effects and the quasi-static 
solutions are not applicable. There is some progress to deal with fully dynamical perturbations by adopting 
a scaling ansatz to solve the 5D equations [64j |. They find that the quasi-static solution is an attractor on 



15 




FIG. 9: The comoving distance r(z) is shonw for ACDM (long dashed), DGP (solid, thick) and the equivalent GR 
dark energy model on the left. On the right, the growth history g(a) — 5(a) /a is shown for LCDM (long dashed) 
and DGP (solid, thick). The growth history for a dark energy model (short dashed) is also shown, with the same 
expansion history as DGP. Due to the time variation of Newton's constant through /3 in Eq. (|53p . the growth factor 
g(a) receives an additional suppression compared with the dark energy model. DGP-4D (solid, thin) shows the 
incorrect result in which the inconsistent assumption is adopted. We set the density parameter for matter today as 
Q m o = 0.3. From [U 




FIG. 10: The galaxy-ISW cross-correlation coefficient R] in each galaxy bin from z — to z — 3. Solid curves denote 
flat ACDM and dashed curves denote open DGP. Note the much larger correlation at high z in open DGP. From [o4| ]. 



subhorizon scales. The ISW effects are sub-dominant compared with the primordial anisotropics formed 
at the last scattering surface. In order to extract the ISW effects, it is proposed to take cross correlation 
between the matter distributions and the CMB [gH, [6(| [67| • It was shown that the quasi-static solution is 
valid to calculate the corss correlation for large I where a signal is maximized (64j. The growth function 
g(a) changes at earlier times in the self-accelerating universe than in the ACDM model. This gives a larger 
signal in the cross correlation at high redshift. Thus higher red-shift galaxies can test the predictions in the 
self-accelerating universe with high significance. 

Hence, structure formation tests are essential for breaking the degeneracy with dark energy models in 
GR [H, [H, [TO, [7l|, • The distance-based SN observations draw only upon the background 4D Friedman 
equation (|4T)]) in DGP models, and therefore there are quintessence models in GR that can produce precisely 
the same SN redshifts as DGP. By contrast, structure formation observations require the 5D perturbations 
in DGP, and one cannot find equivalent GR models. This leads to an exciting possibility to find a failure of 
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w o 

FIG. 11: Equations of state found using two different combinations of data sets. Solid contours are for fits to SN 
la and CMB data, while dashed contours are for fits to weak lensing and CMB data. The significant difference 
(inconsistency) between the equations of state found using these two combinations is a signature of the DGP model. 
The inconsistency is an observational detection of the underlying modified gravity DGP model (assumed here to 
generate the data). From [691 ] . 



GR [69| . Suppose that our Universe is described by the DGP model. However, astronomers still try to fit 
the data by dark energy models in GR. For example, they use the parametrization of the equation of state 
of dark energy 

w = wq + w\z. (56) 

Combining SN observations, CMB shift parameter and weak lensing, there appears an inconsistency. This is 
because weak lensing probes the growth of structure and the growth rate in the DGP model cannot be fitted 
by the growth rate in GR models given the same expansion history. Fig. 1111 demonstrates this possibility. 
In order to quantify the difference in the growth rate, it is convenient to parametrize the growth rate as 

M 

g(a) = cx pj^ dlna(Q(a,y - 1)1 . (57) 

In a quintessence model, 7 is well approximated by 

j(w) = 0.55 + 0.05(1 + w(z = 1)). (58) 

In the DGP, 7 is well approximated as 7 = 0.68 [6(J. Recently, several authors tried to estimate how 
accurately we can constrain 7 using weak lensing in future surveys 0, 0, [75j . These results suggest that in 
the future, we will be able to discriminate ACDM and the DGP model from the difference in the growth rate. 
Fig. [12] shows the constraint on 7 for the DGP model assuming the 'bench mark' survey on weak lensing, 
where the mean redshift is zmean = 0.9 and the number of sources per arcmin 2 is d = 35, 50, 75 flE^ . 

However, as we mentioned before, the weak lensing measure requires knowledge of the non-linear power- 
spectrum. In the DGP, this is a subtle problem. The DGP approaches GR on small scales. This is essential 
to evade the tight constraints from the solar system experiments. The non-linear power spectrum would be 
sensitive to this transition from Brans-Dicke linear theory to GR non-linear theory. The analyses so far have 
used the simple mapping formula developed in GR to derive the non-linear power spectrum. This approach 
could be inconsistent. Nevertheless, the conclusion that we will be able to distinguish the difference in 
the growth of structure would be valid and this is a very exciting possibility that we can achieve in future 
observations. 
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FIG. 12: Confidence regions at 68% for the benchmark survey z mea n — 0.9, d = 35 (outer contour) and for d — 50, 75 
(inner contours) for DGP. The dotted line represents the ACDM value. From fn§ . 



3. Non-linear structure formation 

For quasi-static perturbations, it is possible to extend the linear result to non-linear perturbations by 
taking into account partially the non- linear effects of gravity. A key is the so called brane bending mode 
[49l uA f7e| . This mode describes the perturbations of the location of the brane and mediates an additional 
scalar interaction. In the linear regime, this is the scalar mode that makes the theory of BD type. This 
scalar mode becomes non-linear on much larger scales than gravity. In terms of the brane bending mode tp, 
the effective equations on the brane are given by (7o| 

^V 2 $ = ~8ttG 4 5p + -yVV, (59) 
* + $ = <p, (60) 
where the equation of motion for tp is given by 

3f3(t)Kv + 4 [( V W - (ViV^) 2 ] = 8wG 4 5p. (61) 
a a 

Again these equations are derived by properly solving the 5D equations and imposing the regularity condition 
in the bulk. Here we assume gravity is linear, $ <C 1, but we take into account the second order effects 
of <p. Note that the coefficient of the second order terms is given by As we take r c ~ H^ 1 , the second 
order terms can be comparable to the linear term even if gravity remains linear. 

These non-linear equations are difficult to solve in general. If we assume spherical symmetry, the solution 
for (p is given by 

5 = 3 A «' Mr) = ^\-) i \/:' + (-) - i I- <«*) 




where 



and r g is the Schwarzschild radius r g — 2G±M. The solutions for $ and ^ are obtained as 
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For r > r t , we recover the solutions for linear perturbations 
perturbations are given by [59] 




For r < r*. the solutions for metric 



(66) 



* = -^ + ~ 



(67) 



In this region, the corrections to the solution in 4D GR are suppressed, so that Einstein gravity is recovered. 
The radius r* is the Vainstein radius in the cosmological background. 

The conservation of the energy momentum tensor holds as in GR. Then the continuity equation and the 
Euler equation are the same as in GR: 



OS 



1 



0. 



dvj 
dt 



-(v j Vj)vi 



Hvi = — Vi*, 
a 



(68) 
(69) 



where Mj is the velocity perturbation of dark matter. Eqs. (|59|). (f60|) . (|6Tj) ([68| and (|69|) form a closed 
set of equations that has to be solved to address the non-linear structure formation problem in the DGP 
model. In order to see how GR is recovered dynamically, let us consider the evolution of a spherical top-hat 
perturbation S(t,r) of top-hat radius Rt, where p(t,r) = p(t)(l + 5) is the full density distribution and p(t) 
is the background density [59| . The Newtonian potential 'J dominates the geodesic evolution of overdensity. 
Then the evolution equation for the over-density S is given by 



■■ 4 5 2 
a — — - 



31 + 5 



2HS = AwGp5{l + 6) 



2 1/ /T — 



1) 



(70) 



8r 2 c r g 8(1 + fi ro ) 2 



9/3 2 i?f 



9 (i + n*. 



(71) 



In the linear regime, <5 <C 1, e <C 1, we recover the linear evolution of the overdensity, Eq. (|53|) . On the 
other hand for e ^> 1, the right hand side of Eq. (|70|) becomes the same as in GR and the dynamics of the 
non-linear collapse becomes the same as in GR. Fig. [TBI shows the behaviour of 5 in the DGP compared with 
the ACDM model. 

For the non-spherically symmetric case, we need to solve the equations numerically. We should emphasize 
again that the analysis of the non-linear transition of the theory to GR is essential for the prediction of weak 
lensing and this is an outstanding open problem. 



C. Theoretical consistency of the DGP model 



Although the DGP model offers a concrete example for a modified gravity alternative to dark energy, this 
model is not free from problems. In fact, this model demonstrates the difficulties of modifying GR at large 
distances. One of the problems is related to the non-linearity of the scalar mode. The non-linear interactions 
of the scalar mode become important at the Vainstein length r* ~ { r2 ' r g) 1 ^ '■ If we consider a Planck scale 
mass particle, this length is given by A" 1 = (r^M" 1 )- 1 / 3 , which is - 1000 km for r c - Hq 1 . This defines 
the length below which quantum corrections for the scalar mode become important. Thus A c plays the 
same role as the Planck scale in GR. Then below the length A" 1 , the classical theory loses its predictability. 
This is known as the strong coupling problem fn\ . [TH, ■ There have been debates whether this is indeed 
a problem or not [8(|. It is suggested that there exists a consistent choice of counter-terms for which the 
model remains calculable (78j . 

The most serious problem in this model is that there are ghost-like excitations around the self-accelerating 
universe [zElzEIUli^, [H,[83|. In fact the growth rate already manifests this problem. The solution for the 
linearized perturbations is described by a BD theory with BD parameter given by [4^ . [6l| 



w = |C9-l). 



(72) 
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FIG. 13: Numerical solution of spherical collapse. The left panel shows the evolution for a spherical perturbation 
with Si = 3 x 10 -3 at Zi — 1000 for Q. m = 0.3 in DGP gravity and in ACDM. The right panel shows the ratio of the 
solutions once they are both expressed as a function of their linear density contrasts. From (5i| 



For large Hr c , (3 is always negative. In fact, if u> < —3/2, the BD scalar field has the wrong sign for its kinetic 
term and it becomes a ghost. For de Sitter spacetime, the condition oj < —3/2 implies Hr c > 1/2. This is 
exactly the condition that there exists a ghost in the theory. We can understand the extra suppression of 
the growth rate as due to the repulsive force mediated by the ghost. If we avoid the negative norm state 
when quantizing the theory with ghosts, the ghosts have unboundedly negative energy density and lead to 
the absence of a stable vacuum state. In a Lorentz invariant theory this instability is instant as the decay 
rate of the vacuum is infinity. It is suggested that if there is a Lorentz non-invariant cut-off in the theory and 
the cut-off scale is enough low, it is possible to keep the instability at unobservable level [85|. In the DGP 
model, the strong coupling scale A c may serve as the cut-off scale. It is needed to calculate the decay rate of 
the vacuum and to see whether the self-accelerating universe can survive beyond the age of our Universe. It 
is also necessary to check the validity of the linearized analysis [76|, [H, l87l. [8a. [89] | . Several non-perturbative 
solutions indicate that the self-accelera ting universe would be unstable [88|, [89| . Then we are naturally lead 
to ask what does the solution decay to |9dl |. This is still an open question. See [9l| for a review on the issue 
of the ghost in the DGP model. 

Finally, it was pointed out that time-dependent perturbations around the spherically symmetric spacetime 
have a sound speed greater than 1 92]. Again there are debates whether this is a problem or not. One 
subtlety is that this argument is based on the effective theory for the scalar mode and it is not clear this 
effective theory captures the property of the full gravitational perturbations in the model [93| ■ In addition, 
causality should be defined in a 5D spacetime and it is not clear that the super-luminality in the 4D effective 
theory really means the breakdown of causality in the full 5D theory. 

Although it is still not clear whether we should deny the DGP model as a consistent theory due to these 
problems, this certainly demonstrates the difficulty for the large distance modification of gravity to explain 
the late time accelerated expansion of the Universe. It is necessary to seek improved models that can avoid 
these problems. 



V. CONCLUSION 



In this article, we review the attempts to address the cosmological constant problem and the dark energy 
problem in braneworlds. 
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The cosmological constant problem has resisted solution for many years. The conventional approach 
relies on 4D low-energy physics. This was a natural way of attacking the problem as in a conventional KK 
compactification, the size of extra dimensions must be smaller than TeV -1 and, below the TeV scale, our 
Universe can be described by the 4D effective theory. However, the braneworld picture completely changes 
the notion of extra dimensions. The extra dimensions can be large. For 6D spacetime, the size of the extra 
dimensions can be L <~ 10/zm, with the 6D planck scale 10 TeV. Then above the scale L _1 , the Universe is 
described by 6D and the 4D effective theory cannot be used. In fact the energy density for the cosmological 
constant necessary to explain the present accelerated expansion is roughly p\ ~ L~ 4 . Moreover, the way 
the vacuum energy gravitates in our 4D Universe is completely different in the braneworld. Again in a 6D 
spacetime, the vacuum energy on a 4D brane does not curve the 4D spacetime but just changes the geometry 
of the extra dimensions. This leads to the self-tuning idea where the change of the vacuum energy in 4D 
spacetime is compensated by the modification in the geometry of extra dimensions. Although it was shown 
that the simple non-supersymmetric model does not work, it is hoped that the supersymmetric version of 
the model can realize the self-tuning. A close inspection reveals many problems in this approach, but further 
studies are necessary to judge whether the self-tuning idea really works or not. The outstanding problem 
is to know whether the 4D spacetime settles down to a static solution due to the self-tuning mechanism if 
there is a phase transition in the 4D spacetime. This requires a regularization of the branes and the analysis 
of the time dependent dynamics in the 6D spacetime. Based on the hope for the existence of self-tuning, 
the Supersymmetric Large Extra Dimensions (SLED) model is proposed as a framework to address the 
cosmological constant problem and the dark energy simultaneously. The self-tuning mechanism is supposed 
to cure the problem of the large vacuum energy produced by the phase transition in the 4D spacetime. This 
mechanism relies on the supersymmetry in the 6D spacetime, but supersymmetry is inevitably broken on a 
brane at least at TeV scale. This breakdown is mediated to the bulk only gravitationally and creates a weak 
potential for the radion which is the size of the extra dimensions. The potential energy is determined by the 
supersymmetric breaking scale in the bulk and it is argued that if the size of the extra dimensions is 10^m, 
the 6D Planck scale is 10 TeV and the potential energy for the radion has the right amplitude to explain 
the present accelerated expansion of the Universe. The potential depends on the details of the spectrum 
of theory and it remains to be seen whether this proposal can work or not in a concrete realization of the 
models in string theory. 

The late-time accelerated expansion of the Universe is a new problem forced by the discovery made by 
astronomers in 1998. An interesting possibility to explain this is a large distance modification of gravity. 
Again, the braneworld picture plays an essential role. The braneworld model provides a concrete example 
where gravity leaks off the brane and modifies the 4D GR on the brane at large distances. The DGP model 
is the simplest model that realizes this idea. The action for the model is very simple. The 5D spacetime is 
just a Minkowski spacetime described by Einstein gravity. We are living on a 4D brane where 4D gravity 
is assumed to be induced. Despite the simple set-up of the model, gravity in this model is remarkably 
complicated. In fact there exists a solution (the self-accelerating universe) where the accelerated expansion 
of the Universe is realized just by the modification of gravity. We focused on the possibility to distinguish 
this model from dark energy models in GR by combining various observations. This leads to an interesting 
possibility to find a failure of GR at cosmological scales. Although the DGP model is the simplest model 
where we can address many issues from a simple action, the model is not free from problems. In particular, 
it has been shown that there exists a ghost in a self-accelerating universe. It is crucial to study how we can 
avoid the decay of the self-accelerating universe in order for the observational tests of the model to make 
sense. 

The attempts to use higher-dimensional gravity and branes to address the cosmological constant problem 
and dark energy are new but there has been much progress. In this article, we only covered several attempts 
among them. We see that these attempts bring us a completely new way of attacking long-standing and 
tough problems although none of the models is completely successful so far. We hope further developments 
of the models based on these attempts lead to solutions for the long-standing problems. 
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